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►the  bulk  model  uses  the  entrainment  hypothesis  and  mean- 
turbulent-field  modeling  technique  of  Garwood  (1977).  The 
model  atmosphere  used  in  all  experiments  treats  the  planetary 
boundary  layer  in  a bulk  fashion  and  uses  a highly  simplified 
version  of  the  Arakawa  and  Schubert  (1973)  theory  to  para- 
meterize cumulus  convection.  Model  verification  is  made 
against  data  from  the  undisturbed  period  of  BOMEX  Phase  III. 

Near  the  time  of  maximum  solar  heating,  the  diffusive 
model  of  the  oceanic  mixed  layer  maintains  turbulent  mixing 
down  to  a much  greater  depth  than  the  bulk  model.  In  addition, 
the  diffusive  model  mixes  heat  supplied  near  the  surface 
downward  on  some  finite  time  scale  while  the  bulk  model  assumes 
that  heat  supplied  near  the  surface  is  mixed  instantaneously 
throughout  the  mixed  layer.  For  these  reasons,  the  diffusive 
model  gives  a more  accurate  simulation  than  the  bulk  model 
during  the  day  when  the  net  surface  beat  flux  is  changing 
rapidly.  Both  models  give  good  predictions  of  the  observed 
temperature  profiles  near  sunrise. 

Diurnal  variations  of  cumulus  cloud  activity,  low- 
level  temperature,  moisture,  and  wind  speed  in  the  coupled- 
model  are  found  to  be  forced  primarily  by  direct  radiational 
heating  of  the  atmosphere  rather  than  by  direct  radiational 
heating  of  the  ocean  and  its  associated  diurnal  fluctuation 
in  sea  surface  temperature.  The  direct  radiational  beating 
of  the  atmosphere  forces  the  diurnal  fluctuation  iq  cloud 
activity  through  its  effect  on  the  "cloud  layer  forcing  term" 
of  the  Arakawa-Schubert  theory.  This,  in  turn,  produces  a 
diurnal  variation  in  low-level  moisture  and  wind  speed  through 
cloud-related  vertical  eddy  fluxes  of  moisture  and  momentum  at 
the  top  of  the  atmospheric  mixed  layer. 

— VThe  diurnal  fluctuation  of  the  surface  sensible  plus 
latent  heat  flux  is  controlled  by  both  direct  radiational 
heating  of  the  atmosphere  and  direct  radiational  heating  of 
the  ocean.  However,  the  amplitude  of  this  fluctuation  is 
about  two  orders  of  magnitude  smaller  than  the  amplitude  of 
the  surface  solar  flux  which  is  primarily  responsible  for 
driving  the  diurnal  sea  surface  temperature  fluctuation. 
Although  the  sea  surface  temperature  fluctuation  alters  the 
surface  solar  flux  through  its  influence  on  cloud  cover, 
this  alteration  is  extremely  small.  Consequently,  the  feedback 
loop  between  the  ocean  and  the  atmosphere  is  exceedingly  weak 
at  the  diurnal  time  scale. 
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ABSTRACT 


One-dimensional  coupled  air-sea  models  are  used  to 
compare  the  responses  to  diurnal  heating  of  a diffusive 
model  and  a bulk  model  of  the  oceanic  mixed  layer,  and  to 
investigate  the  diurnal  variability  of  the  upper  ocean  and 
lower  atmosphere  in  the  undisturbed  trade  wind  regime.  The 
diffusive  model  of  the  oceanic  mixed  layer  uses  the  Level-2 
turbulence  closure  theory  of  Mellor  and  Yamada  (1974)  while 
the  bulk  model  uses  the  entrainment  hypothesis  and  mean- 
turbulent-field  modeling  technique  of  Garwood  (1977).  The 
model  atmosphere  used  in  all  experiments  treats  the  planetary 
boundary  layer  in  a bulk  fashion  and  uses  a highly  simplified 
version  of  the  Arakawa  and  Schubert  (1973)  theory  to  para- 
meterize cumulus  convection.  Model  verification  is  made 
against  data  from  the  undisturbed  period  of  BOMEX  Phase  III. 

Near  the  time  of  maximum  solar  heating,  the  diffusive 
model  of  the  oceanic  mixed  layer  maintains  turbulent  mixing 
down  to  a much  greater  depth  than  the  bulk  model.  In 
addition,  the  diffusive  model  mixes  heat  supplied  near  the 
surface  downward  on  some  finite  time  scale  while  the  bulk 
model  assumes  that  beat  supplied  near  the  surface  is  mixed 
instantaneously  throughout  the  mixed  layer.  For  these 
reasons,  the  diffusive  model  gives  a more  accurate  simulation 
than  the  bulk  model  during  the  day  when  the  net  surface  heat 
flux  is  changing  rapidly.  Both  models  give  good  predictions 
of  the  observed  temperature  profiles  near  sunrise. 


Diurnal  variations  of  cumulus  cloud  activity,  low- 
level  temperature,  moisture,  and  wind  speed  in  the  coupled- 
model  are  found  to  be  forced  primarily  by  direct  radiational 
heating  of  the  atmosphere  rather  than  by  direct  radiational 
heating  of  the  ocean  and  its  associated  diurnal  fluctuation 
in  sea  surface  temperature.  The  direct  radiational  heating 
of  the  atmosphere  forces  the  diurnal  fluctuation  in  cloud 
activity  through  its  effect  on  the  "cloud  layer  forcing  term" 
of  the  Arakawa-Schubert  theory.  This,  in  turn,  produces  a 
diurnal  variation  in  low-level  moisture  and  wind  speed 
through  cloud-related  vertical  eddy  fluxes  of  moisture  and 
momentum  at  the  top  of  the  atmospheric  mixed  layer. 

The  diurnal  fluctuation  of  the  surface  sensible  plus 
latent  heat  flux  is  controlled  by  both  direct  radiational 
heating  of  the  atmosphere  and  direct  radiational  heating  of 
the  ocean.  However,  the  amplitude  of  this  fluctuation  is 
about  two  orders  of  magnitude  smaller  than  the  amplitude  of 
the  surface  solar  flux  which  is  primarily  responsible  for 
driving  the  diurnal  sea  surface  temperature  fluctuation. 
Although  the  sea  surface  temperature  fluctuation  alters  the 
surface  solar  flux  through  its  influence  on  cloud  cover, 
this  alteration  is  extremely  small.  Consequently,  the 
feedback  loop  between  the  ocean  and  the  atmosphere  is 
exceedingly  weak  at  the  diurnal  time  scale. 


Section  1 
INTRODUCTION 


The  Naval  Oceanographic  Prediction  System  (OPS) 
has  been  proposed  to  meet  the  growing  need  for  a detailed 
knowledge  of  the  time  evolution  of  the  oceanic  environment 
in  which  various  acoustic  detection  systems  operate.  The 
goal  of  OPS  is  to  provide  five-day  predictions  of  the  major 
features  of  the  world  ocean  with  reasonable  accuracy  and 
detail,  and  two-day  predictions  for  smaller  regions  with  higher 
accuracy  and  detail.  When  it  is  completed,  OPS  will  become 
part  of  the  Primary  Environmental  Processing  System  (PEPS) 
at  Fleet  Numerical  Weather  Central. 

At  any  given  time,  a large  fraction  of  the  world 
ocean  is  covered  by  the  undisturbed  (i.e.  fair  weather) 
trade  wind  regime.  On  the  time  scale  associated  with  OPS, 
the  period  of  the  dominant  fluctuations  in  the  thermodynamic 
structure  of  this  system  is  diurnal.  Furthermore,  as 
noted  by  Shontlng  (1964)  and  Urick  (1975)  Chapter  1,  the 
diurnal  change  in  the  stratification  of  the  upper  ocean  can 
have  an  important  effect  on  the  acoustic  detection  of 
submerged  objects.  Thus,  in  this  paper  we  will  model  the 
diurnal  variability  of  the  upper  ocean  and  lower  atmosphere 
in  the  undisturbed  trade  wind  regime. 

This  study  is  additionally  relevant  to  OPS  because 
the  success  of  the  program  depends  to  a great  extent  on  our 
ability  to  deduce  the  atmospheric  forcing  of  the  ocean  from 
satellite  observations.  Consequently,  it  is  important  to 
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develop  an  understanding  of  the  short-time-scale  variability 
of  processes  on  both  sides  of  the  air-sea  interface  that 
control  the  surface  fluxes  of  heat,  moisture,  and  momentum. 
Also,  because  of  the  relatively  large  diurnal  fluctuation 
in  sea  surface  temperature  characteristic  of  the  undisturbed 
trade  wind  regime,  this  region  represents  an  area  in  which 
interactions  between  the  air  and  sea  could  be  important  on 
the  1-5  day  time  scale. 

We  define  the  marine  boundary  layer  to  be  that 
part  of  the  ocean-atmosphere  system  which  is  directly 
affected  by  conditions  at  the  sea  surface.  In  the  undis- 
turbed trade  wind  regime,  the  marine  boundary  layer  extends 
from  the  base  of  the  oceanic  mixed  layer  (at  about  5 - 50  m 
depth)  to  the  top  of  the  trade  wind  inversion  (at  about 
1.5  - 2.5  km  height).  The  well-defined  layered  structure  of 
this  system  is  illustrated  in  Figure  1 and  has  been  discussed 
by  Garstang  and  Betts  (1974). 

The  oceanic  mixed  layer  owes  its  high  degree  of 
vertical  homogeneity  to  turbulence  that  is  generated  by  the 
wind  and  the  intermittent  upward  buoyancy  flux  through 
the  sea  surface.  A dynamically  stable  water  mass  in  which 
the  vertical  eddy  fluxes  are  very  small  exists  below  the 
mixed  layer.  During  periods  of  relatively  strong  wind 
forcing  and/or  strong  surface  cooling,  the  mixed  layer  tends 
to  deepen  because  the  more  dense  water  from  below  is  entrained 
into  the  layer  by  turbulence  at  its  base.  During  periods  of 
relatively  weak  wind  forcing  and/or  strong  surface  heating, 
however,  the  source  of  turbulent  kinetic  energy  may  become 
too  weak  to  maintain  active  entrainment  at  the  mixed  layer 
base,  causing  the  layer  to  retreat  to  a shallower  depth. 


The  first  attempts  at  modeling  the  upper  ocean 
were  with  diffusive  models.  In  this  type  of  approach, 
vertical  turbulent  fluxes  are  parameterized  in  terms  of  the 
mean  fields  using  eddy  diffusion  coefficients  which  may 
themselves  depend  on  the  mean  fields.  An  advantage  of 
this  type  of  model  is  that  no  assumptions  (other  than 
the  eddy  coefficient  assumption)  need  to  be  made  con- 
cerning the  behavior  of  the  mixed  layer. 

The  first  diffusive  model  of  the  upper  ocean 
was  probably  that  of  Ekman  (1905)  although  this  formu- 
lation was  concerned  only  with  current  structure  and 
took  no  account  of  density  structure.  The  two  problems 
of  current  and  density  structure  were  first  combined  by 
Munk  and  Anderson  (1948).  The  forms  of  the  eddy  coeffi- 
cient in  this  model  were  determined  empirically  from 
stratified  turbulence  data.  Mamayev  (1958)  and  Pandolfo 
(1969)  have  also  proposed  empirical  forms  for  the  eddy 
diffusion  coefficients.  More  recently,  models  such  as 
those  of  Vager  and  Zilitinkevich  (1968)  and  Mellor  and 
Yamada  (1974)  have  been  developed  in  light  of  modern  tur- 
bulence closure  theory.  Martin  (1976)  and  Jacobs  (1978) 
have  performed  numerical  experiments  comparing  several 
of  the  proposed  forms  for  the  eddy  coefficients. 

Another  approach  is  to  make  use  of  the  oceanic 
mixed  layer's  near  vertical  homogeneity  and  model  it  as  a 
"slab"  in  which  mean  field  quantities  have  no  depth  depen- 
dence. This  type  of  treatment  is  advantageous  because  it 
eliminates  the  need  to  solve  for  the  interior  eddy  fluxes  as 
only  the  surface  fluxes  and  the  entrainment  fluxes  at  the 
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layer  base  are  required.  Kraus  and  Turner  (1967)  were  the 
first  to  construct  a slab  or  "bulk"  model  of  the  oceanic 
mixed  layer.  They  used  the  turbulent  kinetic  energy  budget 
for  the  mixed  layer  and  the  approximately  decoupled  state  of 
the  equations  for  thermal  and  mechanical  energies  to  determine 
both  the  entrainment  rate  during  deepening  phase  and  the 
mixed  layer  depth  during  the  retreat  phase. 

The  prototype  model  of  Kraus  and  Turner  (1967) 
has  been  refined  by  several  authors.  Most  of  the  refine- 
ments involve  parameterization  of  the  dissipation  of  turbu- 
lent kinetic  energy  which  was  ignored  in  the  original 
model.  Notable  in  this  regard  are  the  efforts  of  Denman 
(1973),  Niiler  (1975),  Elsberry  et  al.  (1976).  and  Kim 
(1976).  Probably  the  most  significant  alteration  of  the 
protoype  model,  however,  has  been  that  of  Garwood  (1977) 
in  which  mean-turbulent-field  modeling  of  the  vertically 
integrated  equations  for  the  individual  turbulent  kinetic 
energy  components  is  employed  to  close  the  problem. 

Both  diffusive  and  bulk  models  of  the  oceanic 
mixed  layer  may  be  incorporated  into  various  components 
of  OPS.  Therefore,  it  is  important  to  evaluate  and  compare 
their  short-time-scale  responses.  Martin  and  Thompson 
(1978)  have  compared  the  performance  of  diffusive  and  bulk 
models  of  the  oceanic  mixed  layer  on  relatively  long  time 
scales. 


The  theory  which  describes  the  atmospheric  surface 
layer  (see  Fig.  1)  was  first  advanced  by  Monin  and  Obukhov 
(1954).  They  recognized  the  fact  that  the  vertical  eddy 


1-5 


fluxes  of  heat,  moisture,  and  momentum  can  be  regarded 
as  quasi-constant  with  height  in  this  layer  which  implies 
that  the  profiles  of  temperature,  water  vapor  mixing  ratio, 
and  wind  speed  will  have  universal  forms.  This  powerful 
similarity  theory  allows  prediction  of  the  surface  fluxes  in 
terms  of  mean  quantities.  Comprehensive  reviews  of  this 
subject  can  be  found  in  Businger  (1975)  and  Busch  (1977). 

Like  the  oceanic  mixed  layer,  the  atmospheric 
mixed  layer  (see  Fig.  1)  owes  its  tendency  toward  vertical 
homogeneity  to  vertical  mixing  associated  with  its  fully 
turbulent  nature.  In  the  typical  trade  wind  situation,  the 
top  of  the  atmospheric  mixed  layer  occurs  just  below  the 
average  cloud  base  level.  It  is  maintained,  near  this  level 
by  an  approximate  balance  between  the  entrainment  velocity 
at  the  top  of  the  mixed  layer  and  the  mean  downward  motion 
there  caused  by  the  large-scale  mass  divergence  and  the 
compensating  subsidence  produced  by  cumulus  clouds.  Models 
of  the  tropical  marine  atmospheric  mixed  layer  have  been 
proposed  by  Betts  (1973),  Arakawa  and  Schubert  (1974),  and 
Ogura  and  Cho  (1974).  In  all  of  these  studies  the  entrain- 
ment fluxes  at  the  top  of  the  mixed  layer  are  parameterized 
in  terms  of  the  surface  fluxes,  and  the  thin  transition 
layer  which  tops  the  mixed  layer  (see  Fig.  1)  is  treated 
as  a discontinuity  in  the  thermodynamic  mean  field  variables. 

Below  the  transition  layer  the  vertical  eddy 
fluxes  are  accomplished  mainly  by  mechanical  mixing  and  dry 
convection  while  in  the  cloud  layer  above  (see  Fig.  1)  they 
are  due  primarily  to  cumulus  convection.  Parameterization 
of  cumulus  convection  is  one  of  the  most  important  and 
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difficult  problems  in  tropical  meteorology  and  has  received 
much  attention.  Early  attempts  at  cumulus  parameterization 
have  been  reviewed  by  Ogura  (1972)  and  Bates  (1972)  but  the 
most  significant  advance  in  this  area  has  probably  been  the 
spectral  model  of  Arakawa  and  Schubert  (1974).  This  model 
included  important  processes  such  as  compensating  subsidence 
in  the  environment,  entrainment  and  detrainment,  and  the 
re-evaporation  of  cloud  droplets.  It  provided  the  first 
closed  and  theoretically  rigorous,  parameterization  of  the 
interaction  of  a cumulus  ensemble  with  the  large-scale 
environment . 

While  condensation  generally  exceeds  evapora- 
tion in  the  cloud  layer,  the  reverse  is  true  of  the  trade 
inversion  layer  (see  Fig.  1).  The  diagnostic  studies  of 
Nitta  (1975)  and  Soong  and  Ogura  (1976)  indicate  that  clouds 
penetrating  into  the  trade  inversion  detrain  liquid  water 
(which  subsequently  evaporates)  and  saturated  air  creating  a 
strong  heat  sink  and  moisture  source  there.  This  tends  to 
balance  the  heating  and  drying  effects  due  to  large-scale 
subsidence  and  thereby  maintain  the  structure  of  the  inver- 
sion. Finally,  we  note  that  in  undisturbed  trade  wind 
conditions  practically  all  of  the  clouds  originating  from 
the  mixed  layer  are  confined  to  below  the  top  of  the  trade 
inversion  (Nitta,  1975).  Hence,  in  this  case,  the  top  of 
the  trade  inversion  can  be  considered  the  top  of  the  marine 
boundary  layer. 

Because  of  its  relatively  high  degree  of  horizon- 
tal uniformity  and  steadiness,  the  undisturbed  trade 
wind  region  lends  itself  to  one-dimensional  modeling. 

Hence,  following  the  precedents  of  Pandolfo  (1971)  and 
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Pandolfo  and  Jacobs  (1972),  we  will  study  this  regime 
with  a one  dimensional,  coupled  air-sea,  numerical  model. 

In  addition,  we  will  compare  the  performance  of  a dif- 
fusive model  of  the  oceanic  mixed  layer  to  that  of  a bulk 
model  in  simulating  the  response  of  the  upper  ocean  to 
diurnal  heating.  Model  verification  will  be  made  against 
data  from  BOMEX  Phase  III  (see  Kuettner  and  Holland,  1969). 

The  diffusive  and  bulk  oceanic  mixed  layer  models 
are  described  in  Section  2.1.  Section  2.2  develops  the 
atmospheric  model  that  is  coupled  to  both  oceanic  models. 
Section  3 describes  the  experiments  performed  and  Section  4 
discusses  the  results.  Section  5 is  a summary. 
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Section  2 

DESCRIPTION  OF  THE  MODELS 


2.1.  Fundamental  Equations  lor  the  Ocean 

The  fundamental  equations  for  conservation  of  beat, 
salt,  and  momentum  in  the  upper  ocean  can  be  written 


3T 

3t 


- j.7T  - 


-3T 

W3z 


(1) 


3S 

It  = 


(2) 


iX  = 

3t 


fk  x - <x  vz)  X - 


*3^  " w X 


(3) 


where  T is  the  temperature,  S the  salinity,  v the  horizontal 
current  velocity  vector  (standard  right-handed  Cartessian 
coordinate  system  with  z positive  upward  and  origin  at  sea 
surface  assumed),  vg  the  horizontal  geostrophic  current 
velocity  vector,  w the  vertical  component  of  current  velocity. 
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Qr  the  net  radiational  heating  rate,  and  f the  Coriolis 
parameter.  Spatial  averages  at  constant  depth  taken  across 
a region  large  enough  to  encompass  a statistically  signifi- 
cant sampling  of  all  unresolved  phenomena  are  denoted  by 
(“)  and  primes  indicate  departures  from  these  averages. 

Gradients  taken  at  constant  depth  are  denoted  by  (7  ) and 

z 

we  assume  tbat  fluctuations  in  the  horizontal  current  are 
uncorrelated  with  T' , S',  and  v ' . Other  symbols  are  defined 
in  the  Appendix  and  all  of  the  notation  is  standard. 

In  addition,  we  have  a linearized  equation  of 

state 


P 


w 


1 - a 


(T  - Tw)  + B (S  - Sw)| 


(4) 


where  p is  density,  Pw  a reference  density,  Tw  a refer- 
ence temperature,  Sw  a reference  salinity,  and  a and  $ 
empirical  constants. 

2.1.1.  The  Diffusive  Model  of  the  Oceanic  Mixed  Layer 

2. 1.1.1.  Basic  Equations 

In  a diffusive  model,  the  eddy  fluxes  are  para- 
meterized according  to 


where,  for  simplicity  we  have  assumed  j£g  = w = 0.  We 
have  also  assumed  that  the  Rossby  number  of  the  horizon- 
tally averaged  flow  is  <<1,  which  allows  neglect  of  the 
lateral  advection  of  momentum  in  (10).  In  addition,  we 
have  modeled  the  radiational  heating  as 


o 1_ 
p c 9z 
w 


eYZ 


(ID 


where  c is  the  specific  heat  of  seawater,  IG  the  sur- 
face flux  of  penetrating  solar  radiation  (positive 
upward)  and  y an  extinction  coefficient  (see  Niiler 
and  Kraus,  1975).  Furthermore,  we  have  modeled  the 
damping  of  inertial  oscillations  due  to  downward 
propagation  of  internal  wave  energy  (see  Pollard  and 
Millard,  1970)  with  a constant  damping  coefficient 
D in  (10).  As  a final  simplification,  we  will  assume 


KS  - kh 


(12) 


in  all  that  follows. 

2. 1.1. 2 Determination  of  the  Eddy  Diffusion 
Coefficients 

Through  systematic  scaling  of  the  equations  for 
the  Reynolds  fluxes  and  turbulent  kinetic  energy,  Mellor  and 
Yamada  (1974)  developed  a set  of  turbulence  closure  models 
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for  planetary  boundary  layers  labeled,  in  order  of  increas- 
ing complexity,  Levels  1 through  4.  They  used  hypotheses 
proposed  by  Kolmogoroff  (1942)  and  Rotta  (1951)  to  model  the 
triple  correlation  and  dissipation  terms  in  these  equations 
and  the  empirical  constants  arising  from  the  theory  were 
determined  from  neutrally  stratified  turbulence  data. 

In  the  Level  2 model,  which  has  been  applied  to  the 
upper  ocean  by  Mellor  and  Durbin  (1975),  (8)  - (10)  are 
recovered  with  the  eddy  diffusion  coefficients  given  by 


KH 


i 


(13) 


(2  E)  'b  am 


(14) 


where  t is  a turbulence  length  scale  and  E is  the  turbulent 
kinetic  energy 


E 


v’2  + w’2  ) 


(15) 


The  quantities  Ajj  and  Am  are  stability  functions  given  by 


(16) 
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ni  " 


n2  Rf 

d-Rf) 


(17) 


where  nj  - n6  are  empirical  constants  and 


The  turbulent  kinetic  energy  E is  calculated 

from 
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(18) 


The  three  terms  on  the  right  of  (18)  represent  shear  pro- 
duction, buoyant  dissipation,  and  viscous  dissipation  of 
turbulent  kinetic  energy  respectively. 

Finally,  an  estimate  of  the  turbulence  length 
scale,  t , is  calculated  from  the  vertical  extent  of  the 
turbulence  field  according  to 


j[  zE1/2  dz 

/ E1/2  dz 


(19) 


which  closes  the  diffusive  model. 
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Note  that  the  mixed  layer  is  here  defined  as  the 
region  where  E > 0.  Below  the  mixed  layer,  weak  vertical 
diffusion  is  retained  by  setting 


For  a detailed  discussion  of  the  scaling  which  leads  to  this 
model,  see  Mellor  and  Yamada  (1974). 


The  Bulk  Model  of  the  Oceanic  Mixed  Layer 


2. 1.2.1  Basic  equations 


Defining  Tm,  Sm,  and  to  be  the  horizon- 
tally- and  vertically-averaged  mixed-layer  values  of  tempera- 
ture, salinity,  and  current  velocity,  (1)  - (3)  can  be  verti- 
cally averaged  to  yield  conservation  equations  for  these  quan 
ties 


(23) 


where  the  bottom  of  the  mixed  layer  Is  taken  to  be  at  z » -h  and 

t 

( )M  = h y ( 


9 As  before,  we  have  assumed  w = vg  » 0,  neglected  the 

horizontal  advectlon  of  momentum,  modeled  the  radlatlonal 
heating  with  (11),  and  treated  the  damping  of  the  mean 
motion  field  due  to  internal  wave  propagation  with  a constant 
• damping  coefficient. 
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A thin  entrainment  zone  of  thickness  5 is  assumed 
to  exist  at  the  base  of  the  mixed  layer  which  separates  the 
fully  turbulent  water  above  from  the  relatively  quiescent 
water  below.  Assuming  that  the  eddy  fluxes  in  the  entrain- 
ment zone  are  due  entirely  to  the  entrainment  of  the  lower 
fluid  into  the  mixed  layer,  (1)  - (3)  can  be  integrated 
across  the  zone  to  yield  the  familiar  Jump  expressions  for 
the  eddy  fluxes  at  z - -h 

[ • (w' t’ )_h  - -ACw e C ■ T,  S,  ^ (24) 


where  we  is  the  entrainment  velocity  and 

- t_h  - t_h_6  t *=  1,  S,  X (25) 

represent  quasi-discont inuous  changes  in  mean  field  quantities 
across  the  entrainment  zone. 


With  the  mean  vertical  motion  at  the  layer  base 
assumed  zero,  the  relationship  between  we  and  layer  depth 
is  simply 


Thus,  as  the  layer  becomes  shallower  all  entrain- 
ment ceases  and  the  mixed  layer  is  assumed  to  be  completely 
decoupled  from  the  region  below. 

2. 1.2. 2 Entrainment  Hypothesis 

To  close  the  bulk  model,  we  adopt  the  entrainment 
hypothesis  and  mean-turbulent-field  modeling  technique  of 
Garwood  (1977)  which  will  be  summarized  below.  The  basic 
assumption  here  is  that  the  critical  parameter  determining 
the  rate  of  entrainment  is  the  ratio  P of  buoyancy  flux  at 
the  layer  base  (which  stabilizes  the  entrainment  zone)  to 
convergence  of  turbulent  kinetic  energy  flux  there  (which 
destabilizes  the  entrainment  zone).  Thus, 


In  analogy  to  (23),  the  entrainment  buoyancy  flux  can 
be  written 


(b'w' ) = -Abw 

-n  e 


(29) 


(b'w')^  *=  - JaAT  - BASjw 


(30) 


where  the  last  expression  follows  from  (4). 

In  addition,  the  convergence  of  turbulent  kine- 
tic energy  in  the  entrainment  zone  is  modeled  as 


ni  w»(®.  + El)  = Iff ) m 

l (2  ^J-h 


(31) 


where  Te  is  a transport  time  scale  taken  to  be 


Te  = a!h(w'2)  1/2 
e 1 M 


(32) 


2-12 


and  aj  is  a proportionality  constant.  Assigning  the 
critical  constant  parameter  P the  value  of  m4  (following 
Garwood's  notation),  (27)  and  (30)  - (32)  can  be  combined 
to  yield  an  entrainment  equation 


dh 

3t 


ax  h a AT  - 0AsJ 


for 


3h 

3t 


> 


0 


(33) 


To  close  the  model,  the  integral  over  mixed 
layer  depth  of  generation,  redistribution,  and  viscous 
dissipation  of  turbulent  kinetic  energy  are  parameterized 
and  a quasi-steady-state  assumption  for  turbulent  kinetic 
energy  invoked.  Then,  (33)  and  the  budget  equations  for  the 
individual  components  of  turbulent  kinetic  energy  can  be 
solved  for  3h/3t  during  the  deepening  phase.  During  the 
retreat  phase,  the  entrainment  rate  is  identically  zero 
which  allows  h to  be  determined  diagnostically  from  the  same 
set  of  equations.  The  details  of  these  calculations  are 
given  by  Garwood  (1977). 


Finally,  we  note  that  the  heat,  salinity,  and 
momentum  fields  below  the  base  of  the  mixed  layer  are 
again  calculated  from  (8)  - (10)  and  (20). 
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2.2  Fundamental  Equations  for  the  Atmosphere 

The  fundamental  equations  for  conservation  of 
heat,  moisture,  momentum,  and  mass  in  the  lower  atmosphere 
can  be  written 


■i  ■ V - 1 If 


Ip  “'s>  + LA  + <iR 


-i  ■ V “!§  - - % 


(34) 


H - jg)  - (V  • vp)j  - : ^ - |p  .V  (36) 


~ + 3“  - 0 

■*  dp 


(37) 
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where  s is  the  dry  static  energy  given  by 


s » cpT  gz. 


q the  water  vapor  mixing  ratio  (mass  of  water  vapor  per 
mass  of  dry  air),  v the  wind  velocity  vector  projected 
on  a surface  of  constant  pressure,  p the  pressure,  «•>  the 
vertical  p-velocity  (w  «dp/dt),  Qr  the  net  radiational 
heating  rate,  and  Qc  the  rate  of  condensation  minus 
evaporation  per  unit  mass  of  air.  Spatial  averages  at 
constant  pressure  taken  across  a region  large  enough  to 
encompass  a statistically  significant  sampling  of  all 
unresolved  phenomena  are  denoted  by  (—)  and  primes  indicate 
departures  from  these  averages.  Gradients  taken  at  constant 
pressure  are  denoted  by  (vp)  and  we  assume  that  fluctuations 
in  v are  uncorrelated  with  s‘,  q',  and  v'.  Other  symbols 
are  defined  in  the  Appendix  and  all  of  the  notation  is 
standard . 


In  addition,  we  have  the  hydrostatic  equation 


(39) 


and  an  approximate  equation  of  state  for  moist  air 


? - RT( 1 + 0.61q) 
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(40) 


where  p is  the  density  of  moist  air  and  R the  gas  constant. 
2.2.1  The  Atmospheric  Surface  Layer  and  the  Surface  Fluxes 

Under  quasi-steady,  convectively  unstable,  horizon- 
tally homogeneous  conditions,  Monin-Obukhov  similarity 
theory  predicts  that  a surface  layer  of  thickness  -L  exists 
in  which  the  vertical  eddy  fluxes  depart  little  from  their 
surface  values,  and  properly  scaled  vertical  derivatives  of 
temperature,  moisture,  and  momentum  are  universal  functions 
of  z/L.  Here,  z is  height  above  the  surface  and  L is  the 
Monin-Obukhov  length,  given  to  good  approximation  by  (see 
Kraus,  1972) 


pau* 


ly.’-g/Jo  + o.61  (Z^q'  ) 


c T 
P » 


(41) 


where  < is  von  Karmen's  constant,  *a  a reference  density 
for  air,  Ta  a reference  temperature  for  air,  and  u*  the 
surface  friction  velocity  in  the  air  given  by 


(42) 
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with  t the  surface  wind  stress  vector  given  by 


<w'v'  ), 


g 


(43) 


Using  the  fact  that  the  surface  layer  is  a quasi-constant 
flux  layer  and  denoting  quantities  evaluated  at  the  sea 
surface  with  the  subscript  o we  have  to  good  approximation 


w ’ v ’ 


(“’s' )o 
(w’q' )c 
(~’)0 

(~)o 


for  p » p >,  pn  (44) 

o a 


where  p&  is  the  pressure  at  the  top  of  the  surface  layer 
which  is  taken  to  be  at  an  altitude  of  |L|  above  the  surface. 
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The  Monin-Obukhov  similarity  theory  gives  rise  to 
the  following  expressions 


3s  _ ("’s’)0  *s<z/L) 

3z 

PaC 


(45) 


3q 

3Z 


(u,q,)p  V2/L> 
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(46) 
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az 


7£  Vz/L> 


(47) 


where  U is  wind  speed  and  the  <$>  functions  are  empirically 
determined,  non-dimensional  similarity  functions.  The 
forms  of  the  <{)  functions  have  been  reviewed  by  Hogstrom 
(1974)  and  it  is  generally  accepted  that  most  surface  layer 
data  are  well  represented  by 
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for  L < 0 


(48) 
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which  are  the  forms  used  here.  Using  (48),  the  profile  rela- 
tions (45)  - (47)  can  be  integrated  from  z0,  (the  surface 
roughness  length,  to  z (see  Paulson,  1970). 
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(49) 


(50) 


(51) 


where  s0,  qG,  and  UD  are  dry  static  energy,  water 
vapor  mixing  ratio,  and  wind  speed  at  zQ,  and 


v„  - <1  - 9?)1/4 


(52) 


At  the  top  of  the  surface  layer,  S',  q,  and  U are  given 
to  excellent  approximatioj , by  their  mixed  layer  values 
(see  Section  2.2.2).  Thus,  setting  z « -L  in  (49)  - (53) 
yields  after  rearrangement 


Ym  * - 15E>1/4 
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(56) 


where  the  subscript  m indicates  values  characteristic  of 
the  atmospheric  mixed  layer. 

Following  Clarke  (1970)  and  others,  we  assume 
that  the  roughness  length  for  sea  surface  is  given  by 


0.032u*2 

2q  - max  ( . «mln> 


(57) 


where  zm^D  = 0.0015  cm.  We  calculate  sD  from 


s 

o 


(58) 


where  T0  is  the  sea  surface  temperature.  For  the  diffu- 
sive model  T0  is  taken  to  be  the  temperature  at  the 
uppermost  gridpoint  in  the  ocean,  while  for  the  bulk 
model  it  is  set  equal  to  the  mixed  layer  temperature. 

Thus,  we  neglect  any  surface  skin  temperature  effects  (see 
Grassl,  1976).  We  set  qQ  equal  to  0.98  times  the  saturation 
value  at  the  sea  surface  to  account  for  the  effect  of 
salinity  on  vapor  pressure  (see  Roll,  1965).  We  set  UQ 
equal  to  the  magnitude  of  the  surface  drift  velocity  which 
is  approximated  by  (Kraus,  1977) 


v 
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where  (\£g)0  is  the  surface  geostrophic  wind  vector  and 
vm  is  the  mean  current  vector  in  the  oceanic  mixed  layer. 


Equations  (41),  and  (52)  - (57)  are  solved  itera- 
tively at  each  time  step.  Then,  the  surface  sensible  and 
latent  heat  fluxes  (positive  upward)  are  given  by 


-(“'s’  ), 


(59) 
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(60) 


where  Lv  is  the  latent  heat  of  evaporation.  The  surface 
wind  stress  vector  is  assumed  to  be  in  the  direction  of  the 
mixed  layer  wind  velocity  which  implies 


2 *m 


= PftU*  |*m 


(61) 


where  here  vm  indicates  the  wind  velocity  vector  in  the 
atmospheric  mixed  layer. 
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Then,  the  total  surface  eddy  fluxes  of  heat,  salinity, 
and  momentum  for  the  ocean  are  given  by  (see  Niiler  and 
Kraus,  1977) 


Pwc(W'T’)o 


= F 


Bo  " 


+ H_ 


+ L Q 
vo 


(62) 


p (w'S’ ) = S (P  - Q ) 

Kwv  /o  o o o' 


(63) 


"w(wT.>o  = "Jo 


(64) 


where  F is  the  total  solar  radiation  flux,  B the  infrared 
radiation  flux,  I the  penetrating  component  of  the  solar 
flux,  and  P the  precipitation  rate.  The  subscript  o indi- 
cates evaluation  of  quantities  at  the  sea  surface  and  all 
fluxes  are  taken  positive  upward.  Note  that  implicit  in 
(64)  is  the  assumption  of  a fully-developed  sea. 

The  surface  radiative  fluxes  are  related  to  cloud 
cover  using  the  formulation  of  Wyrtki  (1965) 


Fo  = Fi  I 1 " °-38(c  + c2)| 


(65) 


4 */2  2 
B0  = eb^T^O.SB  - 2.0q10  ) (1  - 0.56<T) 


+ 4eb*To  <To  ~ T10> 


(66) 


where  Fj^  is  the  surface  solar  flux  in  the  absence  of  clouds,  C the 
fractional  cloud  cover,  e che  infrared  emissivity  of  the  sea 

surface,  b*  the  Stefan-Boltzmann  constant,  qjo  the  water 
vapor  mixing  ratio  at  10  m height,  and  Tjq  the  air  tempera- 
ture at  10  m height.  Finally,  we  follow  the  suggestion 
of  Niiler  and  Kraus  (1975)  and  take 

I = 0.45  F (67) 

o o 
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2.2.2  The  Atmospheric  Mixed  Layer  and  Transition  Layer 

Defining  sm,  qm,  and  vm  to  be  the  horizontally- 
and  vertically-averaged  mixed-layer  values  of  dry  static  energy, 
water  vapor  mixing  ratio,  and  wind  velocity,  (34)  - (36)  can  be 
vertically  averaged  to  yield  conservation  equations  for  these 
quantities 


~ , 
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where  pb  is  the  pressure  at  the  top  of  the  mixed  layer 
and 


) dp 


(71) 


Note  that  the  Rossby  number  of  the  horizontally 
averaged  flow  is  assumed  <<1  which  allows  neglect  of  the 
lateral  advection  of  momentum  in  (70).  In  addition,  follow- 
ing Arakawa  and  Schubert  (1974)  and  others,  we  have  neglected 
any  condensation  or  evaporation  occurring  within  the  mixed 
layer.  We  mention  here  that  the  assumption  of  vertical 
homogeneity  for  momentum  in  the  atmospheric  mixed  layer  of 
the  trade  wind  regime  may  not  always  be  good  because  of 
baroclinic  effects  (Pennell  and  LeMone,  1974).  However,  for 
simplicity,  we  will  overlook  this  detail  in  the  present 
study. 


The  transition  layer  at  the  top  of  the  atmospheric 
mixed  layer  is  in  certain  ways  analogous  to  the  entrainment 
zone  at  the  base  of  the  oceanic  mixed  layer.  Eddies  in  the 
mixed  layer  impinge  on  the  bottom  of  the  transition  layer, 
overshoot  slightly,  and  thus  entrain  air  from  above  down  into 
the  mixed  layer.  Following  many  others,  we  will  model  the 
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transition  layer  as  a discontinuity  in  mean  field  quantities 
at  Pb.  Thus,  in  analogy  to  (9),  the  eddy  fluxes  at  pb 
can  be  expressed  as 


(-•Ob  " ‘A«“e  « " s»  <*»  X (72) 

where  ue  is  the  entrainment  vertical  p-velocity  and 

mm 

AC  " Cb-  ~ Cm  « * s»  V (73) 


Here  the  subscript  b“  indicates  evaluation  of  a quantity 
just  above  the  discontinuity  at  Pb  (see  Fig.  2)  and 
we  have  assumed  that  Qr  is  continuous  across  the  transi- 
tion layer. 


The  usual  method  for  determining  “e  is  to 
assume  that  it  is  proportional  to  the  surface  buoyancy  flux 
(see  Arakawa  and  Schubert,  1974),  but  we  will  employ  a new 
approach  here.  First,  we  write  a prediction  equation  for 
Pb 


3pb  * 

at"  “ “b  + we  <74> 
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where  horizontal  advection  of  mixed  layer  thickness  has 
been  ignored  and  is  the  average  vertical  p-velocity  in 
the  environment  between  the  clouds.  From  mass  continuity 
we  can  write 


“b  " % + 


(75) 


where  u>b  is  the  vertical  p-velocity  at  pfe  due  to  the 
large-scale  mass  convergence  (see  (37)),  o is  the  fraction 
of  the  model  domain  covered  by  active  cumulus  updrafts,  and 
is  the  convective  vertical  p-velocity  at  pb  defined 
by 


(76) 


with  u*  the  vertical  p-velocity  inside  the  i—  cloud  up- 
draft, the  fraction  of  the  model  domain  covered  by  the 
ith  cloud  updraft,  and  U the  total  number  of  cloud  updrafts. 
In  the  undisturbed  trade  wind  regime,  cr  <<  t a»d  <<  |£>b| 

Thus,  from  (74)  and  (75)  we  have  to  very  good  approximation. 


u 


e 


(77) 
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The  large-scale  vertical  p-velocity  £>b  will  be  regarded 
as  an  imposed  large-scale  forcing  function  for  one-dimen- 
sional implementation  of  the  model.  The  quantity  u*  will 
be  determined  from  a simple  parameterization  given  in  the 
next  section. 

Finally,  following  Betts  (1976)  we  set 


■ Pi 


(78) 


where  pL  is  the  pressure  at  the  lifting  condensation  level. 
The  lifting  condensation  level  is  the  level  at  which  a 
parcel  with  the  thermodynamic  properties  of  the  mixed  layer 
becomes  saturated  when  it  is  raised  adiabatically  from  the 
surface.  In  the  tropics  it  generally  occurs  very  near  cloud 
base.  Thus,  our  conceptual  model  is  that  the  height  of  the 
top  of  the  mixed  layer  is  controlled  by  the  height  of  cloud 
base.  This  is  a good  model  for  the  undisturbed  trade  wind 
regime,  but  it  looses  its  validity  in  subtropical  latitudes 
because  cloud  base  often  falls  below  the  top  of  the  mixed 
layer  there,  yielding  a mixed  layer  that  is  topped  by  a 
solid  deck  of  stratocumulus  clouds  (see  Schubert,  1976). 


2.2.3 


The  Cloud  Layer  and  Trade  Inversion  Layer 


Modeling  the  thermodynamics  of  the  cloud  layer 
and  trade  inversion  layer  is  very  difficult  because,  in 
addition  to  vertical  eddy  fluxes,  subgrid-scale  condensation 
and  evaporation  must  also  be  parameterized.  The  theory  of 
Arakawa  and  Schubert  (1974)  provides  a closed  parameteriza- 
tion of  these  eddy  fluxes  and  source-sink  terms  but,  even 
though  their  model  contains  many  simplifying  assumptions,  it 
is  still  extremely  complicated  and  predicts  more  than  we  are 
interested  in  here. 

A simpler  model  is  appropriate  for  the  short  time 
scale  air-sea  interaction  problem  since  the  main  concern  in 
this  case  is  prediction  of  conditions  near  the  surface 
rather  than  the  details  of  the  interaction  between  the 
cumulus  ensemble  and  the  large-scale  environment  which 
occurs  primarily  in  the  cloud  and  trade  inversion  layers. 

One  way  to  simplify  the  problem  is  to  make  some  assumption 
about  the  pressure  dependence  of  s and  q between  the  top  of 
the  transition  layer  and  the  top  of  the  trade  inversion  and 
integrate  (34)  and  (35)  with  respect  to  pressure  in  this 
region.  The  advantage  of  this  is,  of  course,  that  knowledge 
of  the  specific  pressure  dependence  of  uTsT' , u,"q'* , and  Q 

c 

in  the  cloud  layer  and  trade  inversion  is  unnecessary.  An 
approach  similar  to  this  has  been  taken  by  Albrecht  (1977) 
who  assumed  a linear  variation  with  respect  to  pressure  of  s 
and  q in  the  cloud  layer  and  modeled  the  trade  inversion  as 
a discontinuity  in  s and  q. 
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where  the  subscripts  b~,  c,  and  d indicate  evaluation  of 
quantities  at  the  top  of  the  transition  layer,  base  of  the 
trade  inversion,  and  top  of  the  trade  inversion,  respectively 
(see  Fig.  2)  The  quantities  (3s/3p)I  and  Oq/dp)I  are 
the  vertical  derivatives  of  ss  and  q in  the  trade  inversion 
layer  and,  along  with  s^  and  qg,  will  be  specified  con- 
stants. Thus,  in  a way,  we  follow  the  lead  of  Ogura  and  Cho 
(1974)  by  topping  the  model  with  imposed  profiles  for  s and  q 
In  this  case,  however,  the  region  where  the  profiles  are 
imposed  is  the  trade  inversion  and  not  the  cloud  layer,  as  is 
the  case  for  the  model  of  Ogura  and  Cho.  Imposing  a constant 
stratification  in  the  trade  inversion  is  appropriate  for  the 
problem  considered  here  since  experiments  with  the  model  have 
shown  that  reasonably  large  diurnal  variations  of  the  thermo- 
dynamic characteristics  of  the  trade  inversion  have  a negligi 
ble  impact  on  the  surface  fluxes.  Finally,  we  note  that  the 
assumed  linear  variations  of  s (which  implies  nearly  linear 
variations  of  temperature  and  potential  temperature)  and 
q with  respect  to  pressure  are  good  approximations  to  the 
observations  (e.g.  Augstein  et_  al^.  (1974),  Soong  and  Ogura 
(1976),  Esbensen  (1978)). 


Now,  using  (79)  and  (80),  (34)  and  (35)  can  be  inte 
grated  over  pressure  to  yield 


pressure 


I 


Pb»Pb' 


Pa 


» 


« 


water  vapor  mixing  ratio — ►dry  static  energy— ► 


Figure  2. 


Assumed  model  profiles  of  s and  q.  The 
values  of  p.t  s.f  q.,  /3s\  , and/3g\  are 

W I bpj  I 

imposed.  All  other  quantities  labeled  on 
the  figure  are  predicted  by  the  coupled  model. 
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The  rate  of  condensation  minus  evaporation,  Qc>  can  be 
eliminated  between  (81)  and  (82)  to  yield  a prediction  equa- 
tion for  pc  in  terms  of  the  large-scale  forcing  functions 
(v  . Vs,  v . vq,  u,  Qr  provided  that  Ss^/at,  Sq^/dt  , 

(w's ' )b-,(w 'q' )b~  j (m ’s ' )d , and  (w'q')d  can  be  deter- 
mined. This  will  close  the  model. 

Following  Betts  (1976),  prediction  equations  for 
sj^-  and  q^-  can  be  written 
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where  we  have  made  use  of  (79)  and  assumed  that  the  cloud 
updrafts  at  Pb  have  the  mixed  layer  values  of  s and  q.  We 
have  also  assumed  that  cloud  base  occurs  just  above  Pb  and 
neglected  the  evaporation  of  any  falling  raindrops  there. 
Consequently,  there  are  no  condensation  or  evaporation  terms 
In  (83)  and  (84). 

The  last  terms  on  the  right-hand-side  of  (83) 
and  (84)  represent  important  cooling  and  moistening  effects, 
respectively,  due  to  detrainment  of  cloud  mass  at  Pb 
( w / 3p  < o).  The  diagnostic  studies  of  Betts  (1975)  and 
Nitta  (1975)  have  shown  that  the  detrainment  at  pg-  is  large. 
We  will  parameterize  it  according  to 


(85) 


where  kj  is  a dimensionless  constant.  The  BOMEX  Phase  III 
data  analyzed  by  Nitta  (1975)  suggest  kj  = 3. 

The  eddy  fluxes  at  pg-  and  above  are  assumed  to 
be  due  solely  to  the  presence  of  cumulus  clouds  which  begin 
as  updrafts  with  the  thermodynamic  characteristics  of  the 
mixed  layer.  Thus,  again  invoking  the  assumption  that 
o « 1,  the  eddy  fluxes  of  dry  static  energy  and  moisture 
at  pg- can  be  written  (see  Arakawa  and  Schubert,  1974) 
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(«») ' s ' )^~  *=  -u*  As 

(to^*q ' )b-  = -to*  Aq 


(86) 


The  diagnostic  study  of  Nitta  (1975)  indicates 
that  essentially  all  of  the  cumulus  convection  is  confined 
to  below  the  top  of  the  trade  inversion  for  typical  undis- 
turbed conditions.  Hence,  to  good  approximation  we  can 
write. 


(to's')d  = 0 

(JwT')d  = o 


(87) 
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We  now  turn  our  attention  to  the  convective  vertical 
p-velocity  at  the  top  of  the  mixed  layer,  w*.  In  the 
spectral  model  of  Arakawa  and  Schubert  (1974),  in  which  cloud 
types  are  differentiated  from  one  another  by  their  entrain- 
ment rates,  is  given  by 


r 

o 


m^A)  dA 


(88) 


where  m^(A)dA  is  the  cloud-induced  vertical  mass  flux  at 
Pb  due  to  all  the  clouds  in  the  ensemble  with  an  train- 
ment  rate  of  A, and  Amax  is  the  maximum  entrainment  rate 
for  the  ensemble.  To  calculate  mb(A)  they  introduce  the 
"cloud  work  function,"  A(X),  which  is  an  integral  measure  of 
the  buoyancy  force  for  type  A clouds,  and  demonstrate 
that  a cumulus  ensemble  will  follow  a sequence  of  quasi- 
equilibria with  the  large-scale  atmospheric  forcing.  This 


implies 

Amax 

r 

dA(A)  _ 
dt 

J K( A , A ' ) mb(A')dA’  + F( A ) = 0 

o 

(89) 
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which  is  a Fredholm  integral  equation  that  can  be  solved 
for  mb(M  since  the  kernel,  K(X,X’),  and  the  large-scale 
forcing,  F(X),  can  be  calculated  from  the  large-scale  fields 
alone . 


The  kernel  gives  the  rate  at  which  type  X'  clouds  sta- 
bilize the  environment  for  type  X clouds  per  unit  cloud  mass 
flux  through  release  of  the  moist  convective  instability.  The 
dominant  effect  represented  by  K(X,X')  is  adiabatic  warming  of 
the  cloud  layer  due  to  cloud-induced  subsidence  in  the  environ- 
ment. The  large-scale  forcing  gives  the  rate  at  which  large- 
scale  processes  destabilize  the  environment  for  type  X clouds. 
The  dominant  effects  represented  by  F(X)  in  the  undisturbed 
trades  are  radiative  cooling  in  the  cloud  layer  and  heat  and 
moisture  convergence  in  the  atmospheric  mixed  layer. 


For  the  problem  we  consider  in  this  paper,  the  actual 

form  of  the  spectral  distribution  of  the  cumulus  clouds  is 

not  of  primary  concern.  Instead,  we  require  only  w* 

b 

which  is  an  integral  over  all  the  various  cloud  types.  Thus, 
rather  than  determining  the  spectral  distribution  of  the 
cloud  mass  flux  from  the  full  Arakawa-Schubert  theory,  we  will 
use  some  ideas  from  it  to  parameterize  w*  in  a very  simple 
fashion.  In  particular,  we  write  a simplified  version  of  (89) 
as 


K, 
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(90) 
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where  and  are  some  characteristic  values  of 
the  kernal  and  large-scale  forcing  respectively  and  k2 
is  a proportionality  constant. 

We  will  take  to  be  the  dominant  term  in 
the  kernal  for  non-entraining  (i.e.  A « 0)  clouds  which 
implies 


<pb  - pc 


A 

/ *1 

J ap 


(91) 


where,  for  convenience,  the  average  level  of  detrainment 
has  been  taken  to  be  at  pc.  Similarly,  we  write  F ^ as 


— F + F 
*C  + *M 


(92) 


where  Fq  and  F^  are  characteristic  cloud  layer  forcing 
and  mixed  layer  forcing  terms  respectively.  We  take  Fq  to 
be  the  cloud  layer  forcing  term  from  the  original  theory  for 
non-entraining  clouds,  which  yields 


2-40 


? / 


(-V  ‘V  s 
4 p 


s|£  + QD)  *£ 

<>p  p 


(93) 


Ve  take  Fy  to  be 
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which  is  motivated  by  the  t.rm  arising  from  mixed-layer 
convergence  of  moist  static  energy  in  the  original  theory. 
Hence,  using  (79),  (88),  and  (90)  - (94),  the  convective 
vertical  p-velocity  at  p^  can  be  written 
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A simplified  conservation  equation  for  liquid  water 
(i.e.  cloud  droplets)  can  be  written 


aw 

at 


(96) 


where  W is  defined  by 


N 

w = ]£  w.c. 

i=i  1 1 


(97) 


with  W^  the  average  liquid  water  mixing  ratio  (mass  of  liquid 
water  per  mass  of  air)  of  the  ith  cloud,  C*  the  fraction  of  the 
model  domain  covered  by  the  ith  cloud,  and  N the  total  number  of 
clouds. 


The  first  term  on  the  right-hand-side  of  (96)  can  be 
determined  from  (81)  and  (82)  and  the  precipitation  rate  can  be 
parameterized  according  to 


with  k3  a dimensionless  proportionality  constant.  Thus, 

(96)  can  be  used  to  predict  W,  from  which  an  estimate  of  frac- 
tional cloud  cover  can  be  obtained  by 


(99) 


where  is  a characteristic  value  of  liquid  water  mixing 
ratio  for  a typical  cloud  in  the  ensemble.  Note  that  the 
contribution  to  C of  any  clouds  other  than  those  originating 
as  updrafts  from  the  mixed  layer  are  ignored.  Note  also  that 
o may  be  much  less  than  C since  only  a relatively  small  frac- 
tion of  the  existing  cloud  mass  may  be  characterized  by 
active  updrafts. 

Finally,  we  note  that  no  attempt  will  be  made  to 
model  the  momentum  profile  in  the  cloud  layer  and  above. 
Rather,  the  horizontal  wind  velocity  vector  will  be  imposed 
at  P5  as  an  upper  boundary  condition  for  momentum. 


Section  3 

DESCRIPTION  OF  THE  EXPERIMENTS 


Two  separate  one-dimensional  air-sea  models, 
denoted  by  "D"  and  "B"  are  considered.  The  oceanic  compo- 
nent of  D is  the  diffusive  model  of  Section  2.1.1,  while 
for  B it  is  the  bulk  model  of  Section  2.1.2.  The  atmos- 
pheric component  of  both  models  is  described  in  Section  2.2. 

The  model  equations  (8)  - (10),  (20)  - (22), 

(68)  - (70),  (81)  - (84),  and  (96)  are  solved  using  forward 
time  differencing  with  the  diffusive  terms  treated  implic- 
itly. In  B,  a convective  adjustment  scheme  that  conserves 
heat,  salinity,  and  momentum  is  employed  to  deepen  the  mixed 
layer  in  the  event  that  the  change  in  buoyancy  across  the 
entrainment  zone  (Ab)  becomes  negative.  For  both  D and  B,  a 
uniform  grid  with  spacing  1 m extends  from  the  sea  surface 
to  a depth  of  30.5  m.  A vertical  grid  is  not  required  for 
the  model  atmosphere. 

Values  used  for  the  empirical  constants  needed  in 
the  diffusive  oceanic  model  are  those  given  by  Mellor  and 
Yamada  (1974).  Values  used  for  the  empirical  constants 
required  by  the  bulk  oceanic  model  are  those  given  by 
Garwood  (1977)  plus  P3  = 30  (see  Garwood,  1977).  All 
other  constants  can  be  found  either  in  standard  reference 
books  or  Table  1. 
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Table  1 

Model  Constants  and  Boundary  Conditions 


Quantity 

Symbol 

Value 

Damping  coefficient 
for  inertial  oscillations 

D 

7.5  x 10-6  -1 

Coriolis  parameter 

f 

3.77  x 1C"5  -1 

Empirical  proportionality 
constants  for  model  atmosphere 

kl>  k3>  k3 

3,120,0.15 

Average  liquid  water  mixing 
ratio  of  clouds 

wA 

2 gm  kg-1 

Extinction  coefficient  for 
radiation  in  ocean 

Y 

0.09  m -1 

von  Karman's  constant 

K 

0.4 

Temperature  at  bottom  of 
ocean  models 

Tb 

28 . 025°C 

Salinity  at  bottom  of  ocean 

models 

SB 

34.09  o/oo 

Wind  velocity  vector  at  top 
transition  layer 

of 

(v)b 

-10  i m s-1 

Dry  static  energy  at  top  of 
trade  inversion 

Sd 

3.040  x 10“5j  kg 

Water  vapor  mixing  ratio  at 
of  trade  inversion 

top 

Qd 

5 gm  kg-1 

I 
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We  simulate  Phase  III  Period  A of  BOMEX  (21  - 26  June 
1969)  by  adopting  values  of  the  imposed  large-scale  forcing 
functions  that  are  based  on  observations  during  this  period. 
Table  2 summarizes  the  forcing  functions  for  all  experiments 
reported  here.  Note  that,  following  Pandolfo  and  Jacobs 
(1972),  the  large-scale  vertical  motion  in  the  atmosphere 
and  the  horizontal  advections  will  be  taken  as  time  invariant 
in  the  experiments.  The  forcing  terms  due  to  radiational 
heating  of  the  atmosphere,  however,  are  allowed  to  vary 
diurnally  to  account  for  direct  absorption  of  solar  radia- 
tion by  the  air.  The  amplitude  of  the  heating  rate  is  in 
accordance  with  climatological  estimates  and  the  time-averaged 
heating  rate  is  in  close  agreement  with  the  values  reported 
for  the  BOMEX  area  by  Nitta  and  Esbensen  (1974).  Note  also 
that  the  time-averaged  solar  flux  penetrating  the  sea 
surface  for  all  experiments  is  within  1 - 2 ly  day-*  of 
the  value  observed  at  the  research  vessel  Discoverer  (13° 

08’  N,  53°  51 'W)  during  the  period  in  question. 

We  mention  here  that  Nitta  and  Esbensen  (1974)  ob- 
served a diurnal  fluctuation  in  large-scale  vertical  motion 
over  the  BOMEX  area  which  they  attributed  to  the  land  sea 
contrast.  However,  its  contribution  to  the  diurnal  heating 
rate  of  the  atmosphere  below  the  trade  inversion  is  small 
compared  to  the  diurnal  variation  in  radiative  heating  and 
it  may  be  confined  only  to  the  coastal  regions  anyway. 

Therefore,  it  is  neglected  here. 

The  initial  profiles  of  temperature  and  salinity 
in  the  upper  ocean  are  the  same  for  all  experiments  and 
given  by  observations  at  Discover.  The  initial  velocity 
profile  in  the  upper  ocean  for  all  experiments  is  taken  to 
be  independent  of  depth  and  given  by  V = (-0.08  i + 0.21j)m  s_1. 
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Table  2 

IMPOSED  LARGE-SCALE  FORCING  FUNCTIONS 
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The  initial  conditions  for  the  atmosphere  are  the  same  for 
all  experiments  and  are  obtained  by  imposing  the  large-scale 
forcing  functions  in  Table  2,  holding  the  sea  surface 
temperature  fixed  at  its  initial  value,  and  integrating  the 
model  to  a cyclical  steady  state  from  arbitrary  Initial 
conditions.  Thus,  the  initial  atmospheric  state  is  in 
equilibirum  with  the  large-scale  atmospheric  forcing  and  the 
sea  surface  temperature.  The  time  invariant  upper  and  lower 
boundary  conditions  for  the  coupled  model  are  the  same  for 
all  experiments  and  given  in  Table  1. 

A total  of  four  experiments  are  performed  and 
the  model  configuration  for  each  is  summarized  in  Table  3. 
Experiments  I and  II  are  accomplished  with  the  basic  models 
described  in  Section  2.  Experiment  III  is  the  same  as  I 
except  that  the  sea  surface  temperature  supplied  to  the 
model  atmosphere  is  held  constant  at  its  initial  value 
even  though  changes  are  predicted  by  the  model  ocean.  Thus, 
the  system  is  "half-coupled"  since  the  model  atmosphere  is 
not  allowed  to  "see"  any  changes  in  sea  surface  temperature 
predicted  by  the  model  ocean.  Experiment  IV  is  the  same  as 
I except  that  the  radiational  heating  rate  of  the  atmosphere 
Qr  is  held  constant  at  its  time-averaged  value  for  Experi- 
ment I. 


Table  3 

SUMMARY  OF  EXPERIMENTS 


' I 


Experiment 

I 

II 

III 

IV 


Model 

D 

B 

D 

D 


Description 


Basic  model 


Basic  model 


Half-coupled 


Qr  held  constant 
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Section  4 
RESULTS 


4.1  Mean  Thermodynamic  Structure  of  the  Atmosphere  and 

Mean  Surface  Fluxes 

Fig.  3 shows  vertical  profiles  of  temperature 
and  water  vapor  mixing  ratio  predicted  by  the  model  atmosphere 
averaged  over  the  course  of  Experiment  I.  Also  shown  are 
horizontally-averaged  observed  vertical  profiles  of  these 
quantities  over  the  BOMEX  area  during  the  period  in  question. 
In  general,  agreement  between  the  model  predicted  profiles 
and  the  observations  below  the  base  of  the  trade  inversion 
is  good,  except  for  water  vapor  mixing  ratio  in  the  region 
between  950  and  975  mb.  This  is  because  the  observations  do 
not  show  a high  degree  of  vertical  uniformity  for  q”  in  the 
mixed  layer.  Deardorff  (1978)  has  noted  that  differential 
horizontal  advection  of  moisture  ( i . e. 3/3p(^. Vq ) / 0)  is 
quite  capable  of  destroying  the  well-mixed  structure  of  q 
over  the  tropical  ocean,  and  this  seems  to  be  the  case 
here. 


Note,  however,  from  Figure  3(a)  that  the  well-mixed 
assumption  for  dry  static  energy  seems  to  be  good  in  this 
case  since  the  observed  and  predicted  temperature  profiles 
are  in  extremely  good  agreement  in  the  lowest  50  mb. 

The  greater  tendency  toward  vertical  homogeneity  for  dry 
static  energy  has  been  explained  by  Deardorff  (1978)  in 
terms  of  a feedback  between  the  intensity  of  the  turbulence 
and  the  vertical  gradient  of  dry  static  energy. 
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Because  the  model  atmosphere  does  a good  job  of 
predicting  the  observed  mean  vertical  structure  of  q and  s, 
it  is  expected  that  the  mean  surface  fluxes  of  sensible  and 
latent  heat  are  also  well  predicted.  This  expectation  is 
confirmed  by  Table  4 which  summarizes  various  estimates  of 
H0  + LVQ0  in  the  BOMEX  area  during  the  period  of 
interest.  Note  that  the  average  of  the  estimates  based  on 
observations  is  405  ly  day-*,  while  the  value  resulting 
from  Experiment  I is  415  ly  day-1  which  is  good  agree- 
ment. The  value  of  272  ly  day-1  predicted  by  the  coupled- 
model  of  Pandolfo  and  Jacobs  (1972),  however,  appears 
anomalously  small.  This  may  be  because  all  vertical  eddy 
fluxes  were  parameterized  using  eddy  diffusion  coefficients 
in  their  model  and  no  account  was  taken  of  moisture  phase 
changes  associated  with  cumulus  clouds. 
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Table  4 

Estimated  Average  Surface  Sensible  Plus 
Latent  Heat  Flux  in  the  BOMEX  Area 
for  22  - 26  June  1969 


i (ly  day  _1) 

Source 

Method 

371 

Delmore  (1972) 

oceanic  heat  budget 

379 

Holland  and 
Rasmusson  (1973) 

atmospheric  heat  and 
moisture  budget 

420 

Nitta  and 

Esbensen  (1974) 

atmospheric  heat  and 
moisture  budget 

452 

Nitta  and 

Esbensen  (1974) 

bulk  aerodynamic 
calculations 

272 

Pandolfo  and 
Jacobs  (1972) 

coupled  air-sea  model 
integration 

415 

present  study 

coupled  air-sea  model 
integration  (time  average 
from  Experiment  I) 

4.2  Diurnal  Variability  of  the  Upper  Ocean  and 

Comparison  of  Diffusive  and  Bulk  Models  of  the 
Oceanic  Mixed  Layer 


We  will  now  compare  the  diurnal  response  of  the 
diffusive  and  bulk  models  of  the  upper  ocean.  Fig.  4 
shows  time  series  of  sea  surface  temperature  predicted  by 
Experiments  I and  II  along  with  observations  at  Discoverer. 
Note  here  that,  while  the  surface  fluxes  are  not  exactly 
the  same  for  Experiments  I and  II  because  of  the  difference 
in  predicted  sea  surface  temperature,  at  no  time  does  the 
net  surface  heat  flux  differ  by  more  than  1%  between  the 
two  experiments.  Therefore,  differences  between  the  predicted 
oceanic  fields  of  Experiments  I and  II  are  due  almost 
entirely  to  intrinsic  differences  between  the  diffusive  and 
bulk  models. 

Fig.  4 indicates  that  the  phase  of  the  predicted 
diurnal  fluctuation  in  sea  surface  temperature  is  the 
same  for  both  experiments,  but  the  amplitude  is  larger  for 
Experiment  II.  Furthermore,  good  agreement  is  found  between 
the  results  of  both  experiments  and  the  observations, 
although  the  amplitude  of  the  fluctuation  for  Experiment  I 
seems  to  be  more  realistic  than  that  of  Experiment  II.  Note 
also  from  Fig.  4 that  the  differences  between  the  two  model- 
predicted  curves  are  largest  in  the  afternoon  when  the 
sea  surface  temperature  peaks  but  are  very  small  near  dawn 
when  the  sea  surface  .temperature  is  at  its  minimum.  To 
explain  this  difference  between  the  response  of  the  two 
models,  we  must  examine  the  predicted  internal  structure  of 
the  upper  ocean. 
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Fig.  5 shows  profiles  of  temperature  in  the  upper 
25  m of  the  ocean  for  both  Experiments  I and  II  at  eight 
selected  times  during  the  day  averaged  over  the  five-day 
experimental  period.  Also  shown  are  temperature  profiles 
observed  at  Discoverer  and  averaged  in  the  same  way  by 
Delnore  (1972).  Note  that  the  horizontal  arrow  on  each 
figure  indicates  the  depth  at  which  the  turbulent  kinetic 
energy  predicted  by  model  D goes  to  zero.  Note  also 
that  the  temperature  increase  with  depth  below  the  mixed 
layer  is  characteristic  of  the  BOMEX  area  and  does  not 
represent  an  unstable  stratification  because  of  an  asso- 
ciated salinity  increase  with  depth  in  this  region. 

Fig.  5(a)  shows  the  profiles  at  0800  LST,  two 
hours  after  sunrise,  and  the  mixed  layer  is  still  near  its 
maximum  depth.  Although  the  results  of  both  models  compare 
quite  favorably  to  the  observations,  model  B seems  to  do  a 
better  job  of  predicting  both  mixed  layer  depth  and  tempera- 
ture. 


Fig.  5(b)  shows  the  profiles  at  1100  LST  and  in 
this  case  both  models  show  some  difficulty  in  simulating  the 
observed  temperature  structure.  Moreover,  at  this  point  in 
the  heating  cycle  fundamental  differences  in  the  two  models 
begin  to  appear.  For  example,  in  model  D turbulent  mixing 
extends  down  to  15  m while  model  B predicts  that  all  turbu- 
lence ceases  below  7.5  m.  In  addition,  heat  supplied  near 
the  surface  is  diffused  downward  with  some  finite  diffusive 
time  scale  in  model  D which  results  in  a temperature 


Figure  5(a)-  Vertical  profiles  of  temperature  in  the  upper  ocean  at 
0800  LST  averaged  over  the  6 day  experimental  period.  Observations 

at  Discoverer  (■mhhb)  1 Experiment  I ( ),  Experiment  II  ( ). 

The  arrow  indicates  the  depth  at  which  turbulent  kinetic  energy  pre- 
dicted by  model  D goes  to  zero. 


decrease  with  depth  near  the  surface.  In  model  B,  however, 
any  heat  added  near  the  surface  is  instantaneously  mixed 
throughout  the  mixed  layer. 

Fig.  5(c)  shows  the  profiles  at  1400  LST  and 
here  model  D does  a much  better  job  of  predicting  the 
thermal  structure  than  model  B.  Specifically,  model  B 
predicts  a region  of  excessively  warm  water  between  the 
surface  and  2.5  m depth  and  a region  of  excessively  cool 
water  between  7 and  15  m depth.  This  occurs  because  turbulent 
mixing  near  the  time  of  maximum  solar  heating  is  confined  to 
the  upper  few  meters  in  model  B so  that  heat  supplied  near 
the  surface  is  absorbed  in  a relatively  thin  layer.  Conse- 
quently, this  thin  layer  becomes  anomalously  warm  (hence  the 
larger  diurnal  amplitude  of  sea  surface  temperature  predicted 
by  model  B on  Fig.  4)  while  the  region  below  becomes  anoma- 
lously cool.  For  model  D,  however,  turbulent  mixing  still 
extends  down  to  a depth  of  15  m,  which  seems  to  be  more 
realistic  in  view  of  the  excellent  agreement  between  the 
observations  and  the  predicted  profile  for  D.  The  rate  at 
which  model  D diffuses  heat  downward  during  the  day  is  also 
in  fair  agreement  with  the  observations  of  the  tropical 
diurnal  thermocline  made  by  Ostopoff  and  Northern  (1974). 

In  general,  the  mixed  layer  turbulent  kinetic 
energy  for  bulk  models  changes  instantaneously  with  the 
net  surface  heat  flux.  For  most  diffusive  models,  however, 
the  net  surface  heat  flux  does  not  affect  the  turbulent 
kinetic  energy  throughout  the  mixed  layer  directly. 

Rather,  it  changes  the  turbulent  kinetic  energy  indirectly 
through  its  effect  on  stratification,  which  is  more 
!•*«  - m m 1 1 y realistic  than  the  bulk  model  assumption. 

♦ ,,  • s.-nmtre , since  the  stratification  is  altered  on  some 
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finite  diffusive  time  scale  in  a diffusive  model,  the 
turbulent  kinetic  energy  will  lag  the  net  surface  heat 
flux.  Thus,  if  the  net  surface  heat  flux  changes  sub- 
stantially on  a time  scale  that  is  not  long  compared  to 
the  diffusive  time  scale  of  the  upper  ocean,  then  the 
response  of  a diffusive  model,  because  of  its  inherent 
adjustment  time,  is  expected  to  be  more  realistic  than 
that  of  a bulk  model.  This  expectation  seems  to  have  been 
confirmed  here. 

Fig.  5(d)  - 5(h)  trace  the  thermal  structure 
through  the  rest  of  the  diurnal  cycle.  Note  that  discrep- 
ancies between  the  observations  and  the  predicted  profiles 
for  both  models  D and  B become  progressively  smaller  with 
time  during  the  night.  In  fact,  the  combined  performance  of 
the  models  seems  to  be  best  near  dawn.  This  is  because  the 
net  surface  heat  flux  varies  little  with  time  during  the 
night  which,  in  contrast  to  the  situation  during  the  day, 
allows  the  upper  ocean  to  approach  an  equilibrium  state  with 
the  surface  forcing. 

4.3  Diurnal  Variability  of  the  Lower  Atmosphere  and 

Air-Sea  Interaction 

We  now  consider  the  way  in  which  the  diurnally 
varying  solar  radiative  flux  forces  a diurnal  response  in 
the  model  atmosphere.  The  solar  flux  affects  the  atmosphere 
through  two  mechanisms:  direct  radiational  heating  of  the 
air  and  direct  radiational  heating  of  the  ocean.  The  second 
mechanism  arises  because  the  sea  surface  temperature  provisoes 
the  lower  boundary  conditions  for  temperature  and  humidity  \ 
to  the  atmosphere.  To  examine  the  relative  importance  of 
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these  forcing  mechanisms  and  to  investigate  the  possibility 
of  feedback  processes  between  the  atmosphere  and  the  ocean 
on  the  diurnal  time  scale,  we  will  consider  time  series  of 
various  quantities  predicted  by  the  coupled-model  on  day  6 
of  Experiments  I,  III,  and  IV  (see  Table  3).  Note  that 
because  the  coupled-model  response  is  nearly  cyclical  for 
all  the  experiments,  the  predicted  fields  on  day  6 are 
characteristic  of  any  day  during  the  experimental  period. 

Fig.  6 shows  time  series  of  sea  surface  tempera- 
ture on  day  6 for  Experiments  I,  III,  and  IV.  The  striking 
characteristic  of  this  figure  is  that  there  is  practically 
no  difference  between  the  three  curves.  Since  in  Experiment 
III  the  sea  surface  temperature  supplied  to  the  model 
atmosphere  is  held  constant  at  its  initial  value,  this 
implies  that  changes  induced  in  the  model  atmosphere  by 

these  short  time  scale  variations  in  sea  surface  temperature 
do  not  feed  back  significantly  to  the  model  ocean.  Further- 
more, since  the  radiational  heating  rate  of  the  atmosphere 
is  held  constant  at  its  time-averaged  value  in  Experiment  IV, 
Fig.  6 implies  that  diurnal  fluctuations  in  direct  radia- 
tional heating  of  the  model  atmosphere  do  not  significantly 
affect  the  model  ocean. 

Fig.  7 shows  time  series  of  the  convective  vertical 

p-velocity  at  the  top  of  the  mixed  layer  ut.  Note  that 
* 

takes  on  a maximum  at  night  and  a minimum  during  the 
day  for  both  Experiments  I and  III,  which  is  well  known 
behavior  for  cumulus  convection  over  the  sea  (Hann 
(1901),  Lavoie  (1963),  Kraus  (1963),  Jacobson 
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(1976),  Gray  and  Jacobson  (1977)).  Experiment  IV,  on  the 

other  hand,  exhibits  a much  smaller  diurnal  variation  for 
* 

ujd  and  is  approximately  180°  out  of  phase  to  the  curves 
of  Experiments  I and  III  (see  Fig.  7).  Thus,  since  the  sea 
surface  temperature  supplied  to  the  model  atmosphere  is 
virtually  the  same  for  Experiments  I and  IV  (see  Fig.  6), 
this  implies  that  the  nighttime  maximum  and  daytime  minimum 
in  cumulus  convection  in  the  model  for  Experiments  I and  II 
is  due  to  radiational  heating  of  the  atmosphere.  During  the 
day,  warming  of  the  cloud  layer  by  direct  absorption  of 
solar  radiation  tends  to  decrease  the  cloud  layer  forcing 
term  (see  Equation  93)  and,  as  a result,  the  cloud  activity 
weakens.  At  night,  infrared  cooling  of  the  cloud  layer  has 
the  opposite  effect  and  the  cumulus  convection  increases. 

This  mechanism  is  similar  to  that  proposed  by  Gray  and 
Jacobson  (1977). 

Note  also  from  Figure  7 that  ui£  for  Experiment  I 
is  larger  during  the  day  and  smaller  at  night  than  for 

Experiment  III.  This  is  because  warmer  sea  surface  temperatures 
during  the  day  tend  to  increase  the  cloud  work  function 
and  enhance  the  cumulus  convection,  while  the  reverse  is  true 
at  night  (recall  that  the  sea  surface  temperature  supplied 
to  the  model  atmosphere  is  held  constant  in  Experiment  III). 
Thus,  the  diurnal  sea  surface  temperature  fluctuation  tends 
to  suppress  the  diurnal  variation  in  cumulus  convection  in 
the  model. 

Fig.  8 shows  time  series  of  fractional  cloud 
coverage  C.  For  both  experiments  I and  III,  the  cloud 
coverage  lags  the  convective  vertical  p-velocity  by  about 
8 hours,  thus  taking  on  a maximum  in  the  morning  and  a 
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minimum  in  the  evening.  This  time  lag  reflects  the  fact 
that  individual  cumulus  clouds  do  not  dissipate  immediately 
when  the  updrafts  in  them  cease;  rather,  they  are  character- 
ized by  some  finite  decay  time.  Furthermore,  the  phase 
exhibited  by  experiments  I and  III  is  in  fairly  good  agreement 
with  that  of  the  diurnal  variation  in  radar  echo  occurrences 
over  the  BOMEX  area  reported  by  Hudlow  (1970). 

4t 

As  is  the  case  for  wb,  C for  Experiment  III 
exhibits  a smaller  diurnal  variation  than  the  other  experi- 
ments and  is  approximately  180°  out  of  phase  to  them 
(see  Fig.  8).  This  again  indicates  the  role  of  direct 
radiational  heating  of  the  atmosphere  in  forcing  the  diurnal 
fluctuation  in  cloud  activity.  Fig.  8 also  shows  that  the 
diurnal  amplitude  of  C is  significantly  greater  for  Experi- 
ment III  than  for  Experiment  I which  again  illustrates  that 
the  diurnal  fluctuation  in  sea  surface  temperature  tends  to 
suppress  the  diurnal  variation  in  cloud  activity. 

Fig.  9 shows  time  series  of  precipitation  rate 
P0.  Comparison  of  Figs.  8 and  9 indicates  that  the 
diurnal  behavior  of  the  precipitation  rate  is  qualitatively 
very  similar  to  that  of  the  fractional  cloud  coverage  for 
all  three  experiments.  This  is  a reasonable  result  and  is 
expected  in  view  of  the  parameterization  used  for  PQ  and  C 
(see  Equation  (98)  and  (99)). 


Fig.  10  shows  time  series  of  potential  tempera- 
ture of  the  atmospheric  mixed  layer  0m  (6m  = sm/cp). 

Both  Experiments  I and  III  show  a minimum  in  the  morning  and 
a maximum  in  the  afternoon.  The  curve  for  Experiment  I, 
however,  has  a slightly  larger  amplitude  than  that  of 
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Figure  9.  Time  series  of  precipitation  rate  on 

day  6 for  Experiments  I ( ■) , 1 1 1 ( ), 

and  H ( •) . 
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Experiment  III  and  lags  it  by  about  1 hour.  Thus,  as 
expected,  the  diurnally  varying  sea  surface  temperature 
tends  to  make  6m  slightly  cooler  in  the  morning  and 
slightly  warmer  in  the  evening.  Note  that  the  amplitude  and 
phase  of  the  curve  for  Experiment  I compares  favorably  with 
the  observations  of  the  diurnal  variation  in  temperature  at 

II  m height  reported  by  Paulson  et  al_.  (1972). 

The  curve  for  Experiment  IV  on  Fig.  10  shows  a 
much  smaller  diurnal  variation  than  those  for  Experiments  I 
and  III  and  is  shifted  in  phase  by  about  90°  relative  to 
them.  Thus  we  confirm  that  direct  radiational  heating  of 
the  atmosphere  is  the  mechanism  primarily  responsible  for 
forcing  the  diurnal  fluctuation  in  0m.  This  is  an  expected 
result. 

Fig.  11  shows  time  series  of  mixed  layer  water  vapor 
mixing  ratio  qm.  Note  that  the  curves  for  Experiments  I and 

III  differ  little  from  one  another  with  both  showing  minima  in 
the  morning  and  maxima  in  the  afternoon.  The  amplitude  and 
phase  of  the  curves  for  Experiments  I and  III  are  in  good 
agreement  with  the  observations  of  the  diurnal  variation  in 
water  vapor  mixing  ratio  at  11  m height  reported  by  Paulson 
BX.  ji!.(1972).  The  curve  for  Experiment  IV,  however,  shows  a 
much  smaller  diurnal  variation  with  a maximum  at  night  and  a 
minimum  during  the  day.  Thus,  the  radiational  heating  of 

the  atmosphere  rather  than  the  radiational  heating  of  the 
ocean  (which  produces  a diurnal  variation  in  q0)  is 
responsible  for  forcing  the  diurnal  fluctuation  of  qm  in 
Experiments  I and  III. 
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This  is  an  unexpected  result  since  radiational 
beating,  of  course,  does  not  appear  in  the  conservation 
equation  for  qm.  Consequently,  an  indirect  mechanism 
linking  Qg  to  qm  must  be  present.  This  mechanism  is 
cumulus  convection  which  is  significantly  influenced  by 
Qr  < see  Fig.  7)  and  closely  linked  to  the  vertical  eddy 
flux  of  moisture  at  the  top  of  the  atmospheric  mixed  layer. 
To  illustrate  the  last  point  we  combine  (72)  and  (77)  to 
write 

(*^)b  = (100) 


£ 

Thus,  since  Aq  < 0 and  % < 0,  cumulus  convection  pro- 

duces a flux  of  moisture  out  of  the  atmospheric  mixed 
layer  at  its  top  ((w"^')*,  < 0)  which  tends  to  reduce  qm. 
Recall  that  for  both  Experiments  I and  III,  achieves  a 

maximum  negative  value  near  midnight  and  a minimum  negative 
value  near  midday.  Consequently,  the  atmospheric  mixed 
layer  becomes  dryer  at  night  as  the  enhanced  cumulus  convec- 
tion pumps  moisture  up  into  the  cloud  layer  at  a relatively 
fast  rate.  During  the  day  when  the  cumulus  convection 
becomes  weaker,  moisture  is  still  transported  out  of  the 
atmospheric  mixed  layer  by  the  clouds,  but  at  a slower  rate. 
This  allows  qm  to  increase  in  response  to  the  eddy  flux 
of  moisture  from  the  surface.  Thus,  direct  radiational 
heating  of  the  atmosphere,  acting  through  the  mechanism  of 
cumulus  convection,  controls  the  diurnal  variation  of  mois- 
ture in  the  atmospheric  mixed  layer. 
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Note  here  that  the  diurnal  fluctuation  in  cumulus 
convection  also  produces  a diurnal  fluctuation  in  the  eddy 
flux  of  dry  static  energy  at  the  top  of  the  atmospheric 
mixed  layer  (w's')],).  However,  the  impact  of  this  on  the 
diurnal  variation  of  sm  (and  hence  0m)  is  negligible  in 
Experiments  I and  III  compared  to  the  direct  effect  of  radia- 
tional  heating  in  the  atmospheric  mixed  layer  (i.e.  (QR)ab 
in  Equation  (68)). 

Fig.  12  shows  time  series  of  wind  speed  at  10 
m height  Uiq.  For  both  Experiments  I and  III,  a diurnal 
variation  is  clearly  present  while  for  Experiment  IV  the 
diurnal  signal  is  much  less  evident.  Thus,  direct  radiational 
heating  of  the  atmosphere  is  responsible  for  forcing  the 
diurnal  fluctuation  in  wind  speed  in  Experiments  I and  III. 

The  radiational  heating  of  the  atmosphere  must  affect  the 
wind  speed  through  some  indirect  mechanism  and,  as  was  the 
case  for  q^,  this  indirect  mechanism  is  cumulus  convection. 

The  cumulus  convection  is  linked  to  the  eddy  flux  of  momentum 
at  the  top  of  the  atmospheric  mixed  layer  in  analogous 
fashion  to  the  moisture  flux  (see  Equation  (100)).  This 
eddy  momentum  flux  tends  to  decrease  the  wind  speed  in  the 
atmospheric  mixed  layer  for  the  choice  of  large-scale 
parameters  used  here.  Thus,  since  the  cumulus  convection 
(and,  consequently,  the  eddy  momentum  flux  associated  with 
it)  is  of  greatest  intensity  at  night  for  Experiments  I and 
III,  a minimum  in  low-level  wind  speed  occurs  in  these 
experiments  near  dawn  (see  Fig.  12).  During  the  day  the 
cumulus  convection  weakens  and,  hence,  the  low-level  wind 
speed  increases,  reaching  a maximum  in  the  late  afternoon. 


So,  the  diurnal  fluctuation  in  low-level  wind  speed,  as  was 
the  case  for  the  diurnal  fluctuation  in  qm,  is  controlled 
primarily  by  direct  radiational  heating  of  the  atmosphere 
through  its  influence  on  cumulus  convection. 

Fig.  13  shows  time  series  of  sea  surface  tempera- 
• ture  minus  air  temperature  at  10  m height.  This  quantity 
exhibits  a broad  maximum  in  the  morning  and  a minimum  in  the 
late  afternoon  which  implies  that  the  amplitude  of  the 
diurnal  fluctuation  in  low-level  air  temperature  is  larger 
than  the  amplitude  of  the  diurnal  fluctuation  in  sea  surface 
temperature.  This  is  consistent  with  several  sets  of 
observations  made  in  the  tropics  including  those  of  Paulson 
et  al.  (1972). 


Note  also  from  Fig.  13  that  the  curves  for  all 
three  experiments  differ  substantially  from  one  another  both 
in  amplitude  and  phase.  This  is  also  the  case  for  the  time 
series  of  qQ  - q^o  shown  in  Fig.  14,  which  appear  very 
similar  to  the  curves  of  Fig.  13.  Thus,  direct  radiation 
heating  of  both  the  atmosphere  and  the  ocean  are  found  to  be 
important  in  forcing  the  diurnal  fluctuations  in  sea-air 
temperature  and  moisture  differences  in  Experiment  I. 

Fig.  15  shows  time  series  of  surface  sensible 
heat  flux  H0.  As  is  the  case  for  Tc  - Tjq,  the  three 
curves  differ  substantially  from  one  another,  indicating 
that  direct  radiational  heating  of  both  the  atmosphere  and 
the  ocean  play  an  important  role  in  forcing  the  diurnal 
fluctuation  in  H0.  Note  also  from  Figure  15  that  H0  for 
Experiment  I achieves  a maximum  near  midday  and  a minimum 
near  sunset  and  has  a diurnal  amplitude  of  about  2 ly 
day-*.  The  phase  of  this  diurnal  fluctuation  in  H0  is 
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consistent  with  that  predicted  by  Pandolfo  and  Jacobs  (1972) 
who  also  modeled  BOMEX  Phase  III  with  a coupled  air-sea 
model,  but  the  amplitude  is  only  about  one  half  that  of  the 
earlier  study. 

Fig.  16  shows  time  series  of  surface  latent 
heat  flux  LVQ0.  As  is  the  case  for  H0>  the  three 
curves  differ  substantially  from  one  another,  thus  indicat- 
ing the  importance  of  direct  radiational  heating  of  both  the 
ocean  and  the  atmosphere  in  forcing  the  diurnal  fluctuation 
in  LVQ0.  For  Experiment  I,  the  surface  latent  heat  flux 
attains  a maximum  in  the  early  afternoon  and  a minimum  near 
sunrise  and  has  a diurnal  amplitude  of  about  9 ly  day*. 

As  is  the  case  for  the  surface  sensible  heat  flux,  the 
phase  of  this  fluctuation  is  in  approximate  agreement  with 
that  predicted  by  Pandolfo  and  Jacobs  (1971)  but  the  ampli- 
tude is  only  about  one  half. 

Note  also  from  Figs.  15  and  16  that  the  Bowen 
Ratio  (H0/LVQ0)  averages  about  0.07  for  Experiment  I. 

Holland  (1972)  reports  that  the  observed  Bowen  Ratio  fell 
within  the  range  0.07  - 0.14  during  the  undisturbed  period 
of  BOMEX.  Thus,  our  predicted  value  appears  to  be  too 
small.  The  same  can  be  said  of  the  results  of  Pandolfo 
and  Jacobs  (1972)  who  find  a Bowen  ratio  of  0.05  in  their 
coupled-model  integration  of  the  undisturbed  period  of 
BOMEX. 


Finally,  from  Figs.  15  and  16  we  deduce  that 
the  amplitude  of  the  surface  sensible  plus  latent  heat  flux 
for  Experiment  I is  approximately  9 ly  day-1  which  is 
about  one  half  that  found  by  Pandolfo  and  Jacobs  (1972). 
Delnore  (1972),  on  the  other  hand,  calculated  the  diurnal 
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variation  of  the  surface  sensible  plus  latent  heat  flux  from 
a heat  budget  for  the  upper  ocean  and  found  the  amplitude  to 
be  about  400  ly  day-1  at  Discoverer  during  the  period  in 
question.  The  amplitude  calculated  by  Delnore  for  this 
quantity  appears  unreasonably  large  since  profile  measure- 
ments by  Paulson  £t  aK  (1972)  in  the  BOMEX  area  during  May 
of  1969  did  not  show  these  large  fluctuations.  Rather, 

Paulson  et  al.  (1972)  found  fluctuations  of  H0  + LVQ0 
with  amplitudes  of  about  15  - 20  ly  day  -1 , although  they 
were  not  purely  diurnal.  Delnore's  (1972)  apparent  overesti- 
mate for  the  diurnal  amplitude  of  the  surface  sensible  plus 
latent  heat  flux  may  have  been  partially  due  to  his  neglect  of 
solar  radiation  penetrating  below  the  base  of  the  layer  for 
which  the  heat  budget  was  calculated  and  his  inexact  know- 
ledge of  the  reflectivity  of  the  sea  surface  at  low  sun 
angles.  In  addition,  as  emphasized  by  Stommel  et  al.  (1969), 
internal  waves  on  the  diurnal  thermocline  and  horizontal 
inhomogeneities  in  the  upper  ocean  are  quite  capable  of 
yielding  anomalous  surface  heat  fluxes  from  this  type  of 
budget  calculation. 


Section  5 

SUMMARY  AND  CONCLUSIONS 


We  have  used  coupled  air-sea  models  to  investi- 
gate diurnal  variability  of  the  upper  ocean  and  lower 
atmosphere  in  the  undisturbed  trade  wind  regime.  The 
oceanic  mixed  layer  is  treated  with  two  different  models:  a 
diffusive  model  which  uses  the  closure  technique  of  Mellor 
and  Yamada  (1974)  to  determine  the  eddy  coefficients  and  a 
bulk  model  which  uses  the  entrainment  hypothesis  and  mean- 
turbulent-field  modeling  technique  of  Garwood  (1977).  A 
very  simple  "bulk"  model  of  the  lower  atmosphere  that 
requires  explicit  parameterization  of  the  vertical  eddy 
fluxes  only  at  the  top  and  bottom  of  the  atmospheric  mixed 
layer  is  coupled  to  both  ocean  models.  The  eddy  fluxes  at 
the  top  of  the  atmospheric  mixed  layer  are  parameterized  by 
using  some  of  the  ideas  from  the  model  of  Arakawa  and 
Schubert  (1974);  the  eddy  fluxes  at  the  bottom  of  the 
atmospheric  mixed  layer  are  parameterized  using  Monin-Obukhov 
similarity  theory.  A simple  parameterization  of  rainfall 
rate  closes  the  model  atmosphere. 

Time-independent  values  for  the  horizontal 
advection  of  all  quantities,  the  vertical  component  of  the 
fluid  motion,  the  geostrophic  component  of  the  fluid 
motion,  and  the  lapse  rates  of  dry  static  energy  and  water 
vapor  mixing  ratio  in  the  trade  inversion  are  derived  from 
observations  during  the  undisturbed  period  of  BOMEX  Phase 
III  (21  - 26  June  1969)  and  imposed  as  large-scale  forcing 
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functions.  In  addition,  the  coupled-model  Is  forced  tempor- 
ally with  a prescribed  diurnally  varying  radlatlonal  heating 
profile  which  is  also  based  on  the  BOMEX  data. 

The  model  atmosphere  predicts  the  observed  time- 
averaged  low-level  temperature  profile  quite  well.  The  pre- 
dicted time-averaged  water  vapor  mixing  ratio  profile  also 
compares  quite  favorably  with  the  observations  except  in  the 
upper  two  thirds  of  the  atmospheric  mixed  layer.  This  dis- 
crepancy occurs  because  the  observations  do  not  show  a high 
degree  of  vertical  uniformity  for  water  vapor  mixing  ratio 
in  the  atmospheric  mixed  layer,  while  complete  uniformity  of 
the  quantltly  is  assumed  there  in  the  model. 

The  time-averaged  surface  sensible  plus  latent  heat 
flux  predicted  by  the  model  is  close  to  the  average  of  several 
values  for  this  quantity  calculated  by  various  authors  from  the 
BOMEX  data.  The  model -pred icted  Bowen  Ratio,  however, 
appears  to  be  too  small. 

Both  the  diffusive  and  bulk  models  of  the  oceanic 
mixed  layer  give  good  predictions  of  the  upper  ocean  thermal 
structure  near  sunrise.  This  can  be  attributed  to  the  fact 
that  the  net  surface  heat  flux  changes  slowly  at  night, 
which  allows  the  upper  ocean  to  approach  an  equilibrium  state 
with  the  surface  forcing.  During  the  day,  however,  when  the 
net  surface  heat  flux  changes  rapidly  with  time  because  of 
the  solar  radiative  flux,  the  diffusive  model  performs 
noticeably  better  than  the  bulk  model.  This  is  partly  due  to 
a fundamental  difference  between  diffusive  and  bulk  models: 
a diffusive  model  mixes  heat  supplied  near  the  surface 
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downward  on  some  finite  time  scale  while  a bulk  model 
assumes  that  heat  supplied  near  the  surface  is  mixed  instan- 
taneously throughout  the  mixed  layer.  In  addition,  the  bulk 
model  used  here  predicts  that  all  turbulence  ceases  below  a 
few  meters  depth  near  the  time  of  maximum  solar  heating. 

The  diffusive  model,  on  the  other  hand,  predicts  that 
turbulent  mixing  always  extends  down  to  a depth  of  at  least 
15  m for  the  cases  considered,  which  seems  to  be  a more 
realistic  result  in  view  of  the  observed  temperature  profiles. 

Diurnal  variations  of  cumulus  cloud  activity, 
low  level  temperature,  moisture,  and  wind  speed  in  the 
coupled  model  are  found  to  be  forced  primarily  by  direct 
radiational  heating  of  the  atmosphere  rather  than  by  direct 
radiational  heating  of  the  ocean  and  its  associated  diurnal 
fluctuation  in  sea  surface  temperature.  The  direct  radia- 
tional heating  of  the  atmosphere  forces  the  diurnal  fluctua- 
tion in  cloud  activity  through  its  effect  on  the  "cloud 
layer  forcing  term"  (see  Arakawa  and  Schubert,  1974).  This, 
in  turn,  produces  a diurnal  variation  in  low-level  moisture 
and  wind  speed  through  cloud-related  vertical  eddy  fluxes  of 
moisture  and  momentum  at  the  top  of  the  atmospheric  mixed 
layer . 

The  diurnal  fluctuation  of  the  surface  sensible 
plus  latent  heat  flux  is  controlled  by  both  direct  radiational 
heating  of  the  atmosphere  and  direct  radiational  heating  of 
the  ocean.  However,  the  amplitude  of  this  fluctuation  is 
about  two  orders  of  magnitude  smaller  than  the  amplitude  of 
the  surface  solar  flux  which  is  primarily  responsible  for 
driving  the  diurnal  sea  surface  temperature  fluctuation. 


5-3 


Although  the  sea  surface  temperature  fluctuation  alters  the 
surface  solar  flux  through  Its  Influence  on  cloud  cover, 
this  alteration  is  extremely  small.  Therefore,  the  feedback 
loop  between  the  ocean  and  the  atmosphere  is  exceedingly 
weak  at  the  diurnal  time  scale. 

The  present  study  could  be  extended  in  several  ways. 
First,  an  attempt  could  be  made  to  parameterize  the  evaporation 
of  spray  just  above  the  sea  surface  which  may  be  an  important 
process  in  regions  where  the  low-level  wind  speed  is  strong 
(Ling  and  Kao  (1976),  Wang  and  Street  (1978)).  Neglect  of  this 
process  may  be  part  of  the  reason  why  the  Bowen  Ratio 
predicted  by  the  coupled-model  is  too  small. 

In  addition,  if  the  vertical  eddy  fluxes  of  dry 
static  energy,  water  vapor  mixing  ratio,  and  liquid  water 
mixing  ratio  at  the  base  of  the  trade  inversion  could  be 
parameterized  along  with  the  vertical  distribution  of  the 
raindrop  generation  rate,  the  lapse  rates  of  temperature  and 
humidity  in  the  trade  inversion  could  be  predicted  rather 
than  imposed  quantities.  Although  this  would  make  the  model 
atmosphere  much  more  empirical,  it  would  also  make  it  more 
general.  In  this  case  it  could  provide  an  economical  yet 
detailed  parameterization  of  the  planetary  boundary  layer 
over  the  tropical  ocean  in  an  atmospheric  general  circula- 
tion model. 


5-4 


.VW*. 


APPENDIX 
LIST  OF  SYMBOLS 


A(X) 


aH.  am 


a 


1 


B, 


b 


* 


b 


C 


c 


D 


E 

F(A) 


F 


M 


F 


C 


Cloud  work  function  for  type  \ clouds 

Stability  functions  for  eddy  coefficients 
in  oceanic  mixed  layer 

Proportionality  constant  in  definition  of 
transport  time  scale  for  oceanic  mixed  layer 

Surface  flux  of  infrared  radiation 

Stefan-Boltzmann  Constant 

Buoyancy 

Cloud  cover 

Specific  heat  for  seawater 

Specific  heat  at  constant  pressure  for  air 

Damping  coefficient  for  inertial  oscillations 

Turbulent  kinetic  energy  in  oceanic  mixed  layer 

Large-scale  forcing  for  type  \ clouds 

Characteristic  large-scale  forcing 
for  cumulus  cloud  ensemble 

Contribution  of  mixed-layer  forcing  to  F^ 
Contribution  of  cloud-layer  forcing  to  F^ 


A- 1 


Fq  Flux  of  solar  radiation  penetrating  sea  surface 

F.  Flux  of  solar  radiation  penetrating  sea  surface 

1 in  the  absence  of  clouds 


f 


Coriolis  parameter 


K 

H 

h 

I 


i.  j,  k 

KH,KS,KM 
K(  A. , X ' ) 


Acceleration  of  gravity 

Surface  sensible  heat  flux 

Depth  of  oceanic  mixed  layer 

Surface  flux  of  the  penetrating  component  of 
solar  radiation 

Unit  vectors  in  the  xt  y,  z directions 

Vertical  eddy  diffusion  coefficients  for 
heat,  salinity,  and  momentum  in  ocean 

Kernel  from  Arakawa-Schubert  cumulus  para- 
meterization theory 

Characteristic  kernel  for  cumulus  cloud  ensemble 

Empirical  proportionality  constants  for  model 
atmosphere 

Monin-Obukhov  length 
Latent  heat  of  evaporation 

Turbulence  length  scale  for  oceanic  mixed  layer 


LST 


Local  Standard  Time 


Cloud  base  mass  flux  for  type  clouds 


« 


mb(X) 

« 

m4  Critical  constant  parameter  governing 

entrainment  in  bulk  model  of  oceanic  mixed  layer 


> • 

i 


Total  number  of  active  updrafts  in  cumulus 
ensemble 


t 


t 


r 


ni(i=l,8)  Empirical  constants  in  diffusive  model  of  the 
1 ’ oceanic  mixed  layer 

N Total  number  of  clouds  in  cumulus  ensemble 

P.P  »Pb  Pressure,  pressure  at  top  of  surface  layer,  top 
a of  mixed  layer 

Pb“.PL>Pc  Pressure  at  top  of  transition  layer,  lifting 
condensation  level,  base  of  trade  inversion 

Pressure  at  top  of  trade  inversion 
Precipitation  rate 
Water  vapor  mixing  ratio 


Pd 

Po 

q 


^m’q10  Water  vapor  mixing  ratio  of  atmospheric 

mixed  layer,  at  10  m height 


Qq  Surface  evaporation  rate 

Qr  Net  radiational  heating  rate 

Rate  of  condensation  minus  evaporation  per 
unit  mass  of  air 


<*d 


Water  vapor  mixing  ration  at  model  top 


15 


R Gas  constant  for  dry  air 


A-3 


s 


m 


m 


w 


sd 


T 

T 


10 


VTw 


m 
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AT,  AS  Change  in  temperature,  salinity  across 

entrainment  zone  at  base  of  oceanic  mixed 
layer 


As,  Aq  Change  in  dry  static  energy,  water  vapor  mixing 
ratio  across  transition  layer  at  top  of 
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Change  in  horizontal  velocity  vector  across 
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